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Our understanding of Bethe Anstaz has improved a lot over the last 75 years. This was clear 
from the many excellent lectures on the conference. But there are still unanswered questions and 
actually this lecture will concentrate on four open problems. Two of the problems are related to 
the correlations functions of the XXX spin chain and the XXZ spin chain one to the entropy of 
subsystems and one to the six vertex model with domain wall boundary conditions. 



I. INTRODUCTION 

Originally introduced in 1931 by Hans Bethe in order to solve the isotropic Heisenberg spin chain, the Bethe Ansatz, 
with its numerous generalizations and refinements has been proven an invaluable tool in the field of exactly solved 
models. The versatility of this method is exemplified by the multitude of physical models and problems in which it 
proved useful. Some of these are: J-function Bose gas, massive Thirring model, Hubbard model, XXX, XXZ and the 
XYZ spin chains, six-vertex and eight-vertex models and the list could go on. 

Even though our understanding of exactly solved models has improved a lot there is still a large number of unsolved 
problems. The purpose of this lecture is to present some of them. 

II. THE HEISENBERG SPIN CHAIN AND NUMBER THEORY 

Our first open problem is related to the model that was the starting point for the Bethe Ansatz. Consider the 
antiferromagnetic spin 1/2 Heisenberg XXX spin chain with the hamiltonian 

oo 

n xxx = tfof+i + + -K+i) (i) 

j=—oo 

where af , erf, of are the Pauli matrices and we will denote a the 2x2 unit matrix 

1 \ v ( -i \ z (I \ o / 1 

1 oj aV = [i o) ° ={o -l) a ={o 1 

As we have said this hamiltonian was diagonalized in 1931 by means of what today we call coordinate Bethe Ansatz 
d|. The unique antiferromagnetic ground state in the thermodynamic limit was investigated by Hulthen in 1938 Q 
and the spectrum of excitations which consists of magnons of spin 1/2 was correctly described for the first time by 
Faddeev and Takhtajan 0]. 
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A. Correlation Functions 

The correlation functions of the XXX spin chain are defined as follows: consider n sequential lattice sites and in 
each site pick up a^ k 6 {o^, , af , a 7 ^ } . The average of this operator with respect to the ground state \GS) 

n n 

(II °?) = {GS\\{^\GS) (2) 

k=l k=l 

is called a correlation function. An example is emptiness formation probability which is defined as 

n 

P(n) = (G5|nP 3 -|G5) (3) 
i=i 

where Pj = (cr| + l)/2 is the projector on the state with the spin up in the j-th lattice site. P(n) is the probability 
of formation of a ferromagnetic string of length n in the antiferromagnetic ground state. 

We will show that there is a strong connection between the correlation functions of the XXX spin chain and Riemann 
zeta function at odd arguments, but first we need to introduce some preliminary notions. 

B. Number Theory and Riemann Zeta Function with Odd Arguments 

The roots of all the polynomials with rational coefficients 

r n x n + r„_ix™ _1 + . . . nx + ro = 0, n € Q 

(where we denote by Q the field of rational numbers ) are called algebraic numbers. The transcendental numbers are 
not roots of any polynomials with rational coefficients (i.e. In 2, w, e). We say that x, y are algebraically dependent if 

n— 1 m— 1 



r nm x n y m + r "bx a y b = o, r i:j e 



a=0 6=0 

and algebraically independent if 

n— 1 m— 1 



a=0 6=0 

Also we say that t\ , . . . , t% are algebraically independent transcendental numbers 

r »i.-.»»*r*2 a % = integers, r ai ...a k & Q 

{a.,} 

The Riemann zeta function is defined as [3| 

= E ^7 > 1 ( 4 ) 

n=l 

and it can also be represented as a product with respect to all the prime numbers p (Euler's product) 

«-)=IIt^t ( 5 ) 



p 

p 



For the characterization of the correlation functions it will also be useful to use the alternating zeta series (the value 
of the polylogarithm at root of unity) 



n>0 



( 1 Y"~ - 1 

C»W = Y,-^~ = - U s(-1) (6) 
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where Li s (x) is the polylogarithm. The connection between the Riemann zeta function and the alternating zeta series 
is 

C(S) = Y^lCa( S ) (7) 

where s ^ 1. Unlike the Riemann zeta function which has a pole at s = 1 the alternating zeta function has a limit 
when s — > 1. 

C„(l) = In 2 (8) 

It is known that at even values of the argument the zeta function can be expressed in terms of powers of it but 
at odd arguments the situation is more complicated. R. Apery proved first that £(3) is irrational [B[ and T. Rivoal 
showed that the zeta function at odd arguments takes an infinite number of irrational values @ (see also @]). In 
fact it was conjectured that the values of the zeta function at odd arguments are all irrational even algebraically 
independent transcendent numbers (Don Zagier Q, P. Cartier [§]). 



C. Quantum Correlations and Number Theory 

In 2001 H.Boos and one of the authors were able to calculate exactly the emptiness formation probability P(n) for 
small strings (n — 1, . . . , 4) obtaining the following results 

P(l) = i = 0.5, 
P(2) = i-iln2 

P(3) = \ -In 2 + |c(3) = 0.007624158, 
4 8 

P(4) = ^21n2 + ^ C (3)-^ln2. C (3)-|c 2 (3) 
55 85 
"24 C(5) + 24 hl2 ' C(5) = °- 000206270 - 

P(l) is obvious from symmetry and P(2) can be obtained from the result of Hulthen [2j for the ground state energy. 
P(3) can be obtained from Takahashi's result [l(| for the nearest neighbor correlation (see also [rTf) and P(4) was 
obtained from from the i nteg ral representation derived in [l2j based on the vertex operator approach |70l] . Subsequent 
computations of P(5) in [15| and P(6) in [14( showed that the emptiness formation probability for small strings share 
the same structure. This led to the following conjecture 

Conjecture. 1 (Boos, Korepin 2001) Any correlation function of the XXX spin chain can be represented as a poly- 
nomial in In 2 and values of Riemann zeta function at odd arguments with rational coefficients fZl /. fZ3 /. 

The conjecture was sustained also by computations of different correlation functions such as two-point spin-spin 
correlators some examples being listed bellow 

= - |c„(l) = -0.147715726853315 

(S z 3 S z J+2 ) = ^ - |C«(1) + C«(3) = 0.060679769956435 
1 74 ^fi R 

( S i S Us) = T^- 3 Ca(l) + yCa(3)-yCa(l)Ca(3)--Ca(3) 2 

-yCa(5) + yCa(l)Ca(5) = -0.050248627257235 

(SjSj-h) = T\-fCa(l) + ^C Q (3)~72C Q (l)C a (3)-^C Q (3) 2 -^C a (5) 
4640 , , ^ , , 220 . \ > / \ 400 > , N? 

+ — Ca(l)Ca(5) - — Ca(3)Ca(5) - — Ca(5) 2 
+ ^Ca(7) - ^°Ca(l)Ca(7) + 280C Q (3)Ca(7) 

= 0.034652776982728 
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Again the nearest neighbor correlator was obtained from Hulthen result [2J and the second-neighbor correlator was 
obtained by Takahashi in 1977 [10| using the strong coupling expansion of the ground state energy of the half-filled 
Hubbard chain. The next nearest correlators were calculated in [16j|.[l7j|. 

The conjecture was finally proved in 2006 by H. Boos, M. Jimbo, T. Miwa, F. Smirnov and Y. Takeyama [l8l ]. 
However we still not have explicit expressions for the rational coefficients that enter in the expression for the correlation 
functions as polynomials in In 2 and Riemann zeta function at odd arguments. Such a description will provide a 
tractable formula for the correlations of the XXX spin chain. 

Open problem: Efficient description of the rational coefficients which appear in the expression for 
the correlation functions of the XXX spin chain as a polynomial in alternating zeta series at odd 
arguments. 

III. ENTROPY OF SUBSYSTEMS 

We are interested in the following physical situation. Consider a one-dimensional system of interacting spins (we 
can also consider particles interacting) in the ground-state denoted by \GS). We will treat the whole ground state as 
a binary system A and B where A is a block of neighboring spins and B is the rest of the spins in the ground state. 
The density matrix of the entire system is 

Paub = \GS){GS\ (9) 

and the density matrix of the subsystem A obtained by tracing away the B degrees of freedom is 

PA = Tr B (pA&cB) (10) 

In a seminal paper fl9j Bennet, Bernstein, Popescu and Schumacher discovered that the von Neumann entropy of a 
subsystem A 

S(p A ) = -Tr(p A lnp A ) (11) 

is a measure of entanglement. Entanglement is the fundamental resource used in quantum computation and quantum 
information. A better understanding of entanglement will provide further insight in the theory of quantum phase 
transitions but also in the physics of strongly correlated quantum systems where deeply entangled ground-states play 
a major role in the understanding of these quantum collective phenomena. Consequently a large amount of effort was 
invested in studying the entropy of subsystems in a large class of quantum systems. 

A. General results 

If we consider the doubling scaling limit in which the size of the block of spins is much larger than one but much 
smaller than the length of the entire chain we can present some general results about the behavior of the entropy of 
subsystems. 

In the case of ID critical models (gap-less) the entropy of the subsystems scales logarithmically with the size of the 
block. More precisely for a block of n spins we have 

c 

S(n) = — Inn — > oo (12) 

where c is the central charge of the associated conformal field theory that describes the critical model. This formula 
was first derived for the geometrical entropy (the analogous of pip for conformal field theory) by Holzhey, Larsen 
and Wilczek in [2(| (see also [21] , [22] , [28] ) . Some examples are: Hubbard model , XXO (or isotropic XY) spin chain, 
higher spin generalization of the isotropic XXX antiferromagnetic spin chain, Bose gas with S interaction etc. 

In the case of non-critical models (gap-full) it was conjectured in [2l| (based on numerical evidence for some spin 
chains) that the entropy of subsystems will increase with the size of the subsytem until it will reach a limiting value 
S(oo). This was proved for the XY spin chain in [30j,[3l| where the limiting value was analytical computed and 
also checked for the spin chain introduced by Affleck, Kennedy, Lieb, Tasaki ([23| AKLT model) in [24[ where it was 
showed that S(oo) = 2. This explains why the Density Matrix Renormalization Group (25[(DMRG) technique works 
so well for the AKLT model and fails to reproduce quantum critical behavior. It was noted in [2lT ] that in order for 
the DMRG technique to work we need to have for the S(p A ) a limiting value as the size of the subsystem increases 
(this is equivalent with a bounded rank for p A ) and as we have shown this is not the case in critical models. 
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B. An example: the XY spin chain 



We can make all these general considerations more exact by presenting the results for the XY spin chain in magnetic 
field with the hamiltonian 



7~Lxy 



■ £ ((l + 7)ofof+i + (1- 7)4^1 +k: 

j=~OG 



where < 7 < 1 is the anisotropy parameter an h is the magnetic field. The model was solved by E.H. Licb, T. Schulz 
and D. Mattis in zero magnetic field case 32j and by E. Barouch and B.M. McCoy in the presence of a constant 
magnetic field [26]. The ground state is unique and in general there is a gap in the spectrum. 
The density matrix of a block of n neighboring spins in the ground state can be expressed as 



p(n) 



1 f n \ " 

Yn E [li-T )<GS\I[o?\GS> 

3=1.. n 



and the limiting value of the entropy in the double scaling limit depends on the isotropy and magnetic field. We can 
distinguish three cases: 



• Case la: moderate magnetic field 2^/1 — 7 2 < h < 2 

• Case lb: weak magnetic field including zero magnetic field < h < 2^J 1 — 7 2 

• Case II: strong magnetic field h > 2 

The result for these regions obtained in [30l | (we denote by S(oo) the limiting value of S(p(n)) when n —> 00) is 



/>OC 

5(oo) = - jf In 



) 3 (ix+^)9 3 (ix-^) 



dx 



J z\ — ) / sinh (nx) 

where the modulus k of the theta function is different in the three regions; r = I(k')/I(k) where I(k) is the complete 
elliptic integral of modulus k. k! = yl — k 2 is the complementary modulus and a — 1 in Case I and a = in Case 
II. Using the approach of [28| I. Peschel independently calculated and simplified the results in region (la) and (II) 
obtaining 



5(oo) 



In 



/ k 2 



\16k' 



1 - 



U{k)I{k') 



In 2 



ith 



k = 



v/(V2) 2 +7 2 -l 



(la) 



5(oo) 



12 



with 



(II) 



and the simplified result for the region (lb) was finally obtained in [3C 



5(oo) = i 



in 



/ k 2 



V 16fc' 



i ]f\ 4I(k)I(k' ') 



In 2 with k 



yT - 7 2 - (h/2f 
yfl (h/2) 2 



(lb) 



In the isotropic case 7 = the ground-state is again unique but now the model is critical for h < 2. We expect 
that the entropy will scale logarithmically and indeed it was showed in [27j that 



at 1 ^ l , f n — t jA e 1 cosht/2 1 

S(p(n)) = - ln(nV4 - h 2 ) - / dt i 1 5 ^ } 

VFV " 3 1 Wo I 3t tsinh 2 (t/2) 2sinh 3 (t/2)J 



We see that the result of applying the magnetic field is very simple effectively reducing the size of the subsystem. If 
the magnetic field is larger than the critical value 2 then the ground-state is ferromagnetic and the entropy is zero. 
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C. The XXZ spin chain 



It will be highly desirable if we would be able to obtain the same amount of information about the XXZ spin chain 
with the hamiltonian 



Hxxz = - (^+i+ CT M+i + Ao *°i+i) 

j=—oo 

If A > 1 the ground-state is ferromagnetic so S(p(n)) = and in the critical region (—1 < A < 1) the entropy will 
scale logarithmically (see fT2|) . In the gapped antiferromagnetic case (A < —1) we expect that S(p(n)) will tend to a 
limiting value S(oo) but at this moment the analytic evaluation of this constant is missing. 

Open problem: For the XXZ spin chain in the antiferromagnetic region (A < — 1) calculate the 
limiting value of the entropy of a subsystem when the number of spins in the block is large. 



IV. 



ASYMPTOTICS OF TIME AND TEMPERATURE DEPENDENT CORRELATION FUNCTIONS 



The third open problem that are we going to present is related to the asymptotic behavior of correlations functions 
when the space and time separation is large. As in the previous section we are going to present two examples for 
which we have satisfactory results and which we believe are going to facilitate the understanding of the problem at 
hand. 

Consider the isotropic XY model [32] in transverse magnetic field with the hamiltonian 



where a are Pauli matrices and h is the magnetic field. We are interested in the asymptotic behavior of the time and 
temperature correlation function 



g(n,t) 



Tr {( 



,-n XY /T\ n + 



<(*2)<r-(ti)} 



Tr ( e -*W r ) 



(13) 



when n = ri2 — n\ and t = ti — t\ are large and h € [0,2). In [35[ it was showed that g(n,t) decays exponentially 
but the rate of decay depends on the direction cj) defined as cot0 = n/At when n,t — > oo. The asymptotics in the 
space-like and time-like regions are: 

• Space-like directions < <fi < ir/4 



g(n, t) — ► C exp 



2tt 



dp In 



tanh 



h — 2 cosp 
T 



(14) 



Time- like directions 7r/4 < <j> < ir/2 



l(n,t) -> Ct {2y + +2u - ] exp \ — / dp \n - 4isinp| In 



tanh 



h — 2 cos p 
f 



(15) 



with 



tanh 



h — 2 cos po 
f 



1 

2^ 



tanh 



h + 2 cos po 
f 



— = sin On 

it F 



At zero temperature the asymptotics of the correlation functions were evaluated in [33j | and [34j and for 4> = n/2 
the leading factor was computed in [37j |. 

Similar formulae were also obtained for the ^-function Bose gas (N bosons interacting via a repulsive ^-function 
potential of strength c) which is characterized by the hamiltonian 
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d 2 

^ = -E&2+ 2c E 8(x 3 --x k ) (16) 

3 = 1 3 N>j>k>l 

In this case the asymptotic formula (x, i — ► oo) obtained in [36j using the determinant approach to quantum correlation 
functions lL0| is 



1 /"°° dX 

(i>(0,0)^(x,t)) T -^exp<j — y |:c - u(A)t| In 



tanh 



e(A) 



(17) 



2T 

where e(X),pt{X) and w(A) are the finite temperature dressed energy, pseudoparticle density and Fermi velocity. 

A. The XXZ Spin Chain 

In the last years remarkable progress was made in obtaining multiple integral representation for the correlation 
functions of the XXZ spin chain. We remind the hamiltonian 

oo 

n xxz = - (^+i + *K+i + Ao i°i+i + K) 

J=-O0 

where A is the anisotropy and h the magnetic field. U sing the q- vertex operator approach and corner transfer matrices, 
Jimbo, Miki, Miwa and Nakayashiki obtained in 1992 [38|] multiple integral representation for the correlation functions 
in the massive regime (A < —1), and a conjecture for the critical regime (|A| < 1) was proposed in [3!|. The next 
important steps were made by Kitanine, Maillet, Terras and later Slavnov when they proved the previous results and 
their extension in the case of a magnetic field in both regimes in 1999 [lO, El|. Their method used the algebraic 
Bethe Ansatz [7l| ■ Later they were able to obtain representations for the two-point correlations in terms of a multiple 
integral [42j which they called master equation. The generalization of the method for the time-dependent correlations 
was made in [43[ (see [44| for a review of the entire series of papers) . The final step was the extension of the multiple 
integral representation for finite temperatures which was done by Gohmann, Kliimper, Seel and Hasenclever in [45j,[46[. 
It should also be mentioned that in the study of the correlation functions in integrable models an important role is 
played by the quantum Knizhnik-Zamolodchikov equation [48| (see also [49| . [69T | ) a nd for supersymmetric fermion 
models determinant representation of the correlations functions were reviewed in [471 ] . 

Despite these considerable efforts the asymptotics of time and temperature dependent correlation functions, such 
as g(n, t) (see IT3|) for the XXZ spin chain, are still out of our reach. We expect that they will decay exponentially 
like in the previous examples of the isotropic XY spin chain or the <5-function Bose gas but at this moment we do not 
have such a result. 

Open problem: Derive an explicit analytic formula describing the asymptotic exponential decay of 
time and temperature correlation functions in the critical region (1 < A < 1) of the XXZ spin chain. 



V. SIX VERTEX MODEL WITH DOMAIN WALL BOUNDARY CONDITIONS 

The six vertex model with domain wall boundary conditions (DWBC) was introduced in [50| where recursion 
relations for the partition function were also derived. By domain wall boundary conditions we understand the 
situation in which the arrows on the upper and lower boundaries point in the square lattice and the ones on the left 
and right boundaries point out. 

Using the recursion relations Izergin was able to derive a determinantal formula for the partition function. In the 
case of the homogeneous system if we parametrize the weights as 



a = sin(7 — t) b = sin( , y + t) c—sin(2j), \t\ < 7 (18) 
then the result obtained in [5 11 ] is 

[smb +i)»n(7-l)]"' . , 

Zn = — — * { ' 
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where tn is the Hankel determinant 



and 



ji+fc-2 

^=det(— (20) 

Ub 1 Ki.k<N 



m = . sir f ^ (21) 

v ; sin( 7 + 1) sin( 7 -t) v ' 



A. Connection to Algebraic Combinatorics 



One important feature of the six vertex model with DWBC on a N by N lattice is the connection with the 
alternating sign matrices of the same dimension. An alternating sign matrix (ASM) is a matrix in which all the 
entries are —1,0,1, the sequences of l's and — l's are alternating and the sum of all the elements on each row and 
column is 1. An example of dimension 4 is 



/ 


1 





°\ 








1 





1 


-1 





1 


V o 


1 





0/ 



(22) 



The study of ASM is the domain of algebraic combinatorics, 
in H Iff to be 



The number of ASM of dimension N was conjectured 



N-l 



a(n) = n 



(3n+l)\n\ 
(2n)!(2n + l)! 



(23) 



and it was proved for the first time by Zeilberger [57|. Using the fact that there is a one-to-one connection between 
the allowed configurations of the six vertex with DWBC and ASM and the determinantal formula (fl"9|) Kuperberg 
was able to give an alternative proof 56] to the ASM conjecture (see also [HI]). An unified and simplified treatment 
of ASM enumerations and the relation with the classical orthogonal polynomials was found by Colomo and Pronko 
in (Hi. 



B. The Thermodynamic Limit 



From the physicist's point of view a more important feature of the six vertex model with DWBC is the dependence 
of the bulk free energy on the boundary conditions. The fact that particularly boundary conditions affect the bulk free 
energy of the six vertex model was discovered in [52fl. The study of the thermodynamic limit of the bulk free energy for 
the six vertex model with DWBC was initiated in [531 ] . Using the determinant representation of the partition function 
a Toda differential equation was asymptotically solved in order to extract the bulk free energy in the disordered and 
the ferroelectric phases. It was noticed that for the disordered region the result is different from the result derived 
for periodic boundary conditions but the phase transitions take place in the same place. In (59| P. Zinn- Justin was 
able to rewrite the determinantal formula in terms of a partition function of a random matrix model. This allowed 
him to obtain also the bulk free energy in the antiferroelectric phase and to conjecture that in the disordered region 
the partition function behaves like 



Z N ~ CN K e N2f as N -> oo (24) 
with C a constant and the bulk free energy is 
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Usin g th e results from , Bleher and Fokin obtained the large iV asymptotics of the free energy in the disordered 
region [601 ] . They have used completely integrable integral operators, the Ricmann-Hilbcrt approach and the Dcift- 
Zhou nonlinear steepest descent method. This allowed Bleher and Fokin to prove Zinn-Justin's conjecture obtaining 
for the k exponent (see |2~41 the value 

K = 12 " Mir - 2 7 ) (26) 
C. Correlation Functions 

As in all exactly solvable models the study of correlation functions is of extreme importance. In the case of the 
six vertex model with DWBC an added difficulty is constituted by the lack of translation invariance. However in the 
vicinity of the boundary some results are known. Determinant representations for the 1-point boundary correlation 
function were derived in [65[ using the algebraic Bethe Ansatz. This result was rederived using a direct combinatorial 
approach and extended to 2-point boundary correlation functions in [6l|. In principle this technique can be used to 
express more general boundary correlation functions as sums over determinants. For a different calculation of the 
2-point boundary correlation function and the relation with the doubly refined x-enumerations of ASM see |64| . 

Open problem: Calculate the bulk and boundary correlation functions of the six vertex model with 
domain wall boundary conditions. 

D. Further Generalizations 

The main feature of formula (fl9|) is the fact that the partition function of the six vertex model with DWBC is 
expressed as a determinant of the size of the system. The next obvious step would be the generalization of DWBC in 
the case of other models such that the partition function share the same feature. In the case of the six vertex model 
the generalization for spin-fc/2 was done in [(32j using a fusion procedure, for the level-1 so(n) affine vertex models 
the generalization was made in [o3| and recently DWBC were obtained for supersymmetric models (66j. 

Open problem: Obtain generalizations of domain wall boundary conditions for other models such 
that the partition function is expressed as a determinant of a dimension proportional with the size of 
the system. 

VI. CONCLUSIONS 

Of course the open problems that we have presented do not exhaust the list of unanswered questions in the field of 
exactly solvable models. We should mention for example the solving of the time-dependent Schrodinger equation for 
spin chains which plays an important role in quantum information theory [o7| . In the light of the talks given at the 
conference we can definitely say that Bethe Ansatz is alive and well after 75 years of developments and will carry on 
into the 21st century. 
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